Transposable elements regulate thymus development and function.

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/ß. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

Autophagy degrades immunogenic endogenous retroelements induced by 5-azacytidine in acute myeloid leukemia.

The hypomethylating agent 5-azacytidine (AZA) is the first-line treatment for AML patients unfit for intensive chemotherapy. The effect of AZA results in part from T-cell cytotoxic responses against MHC-I-associated peptides (MAPs) deriving from hypermethylated genomic regions such as cancer-testis antigens (CTAs), or endogenous retroelements (EREs). However, evidence supporting higher ERE MAPs presentation after AZA treatment is lacking. Therefore, using proteogenomics, we examined the impact of AZA on the repertoire of MAPs and their source transcripts. AZA-treated AML upregulated both CTA and ERE transcripts, but only CTA MAPs were presented at greater levels. Upregulated ERE transcripts triggered innate immune responses against double-stranded RNAs but were degraded by autophagy, and not processed into MAPs. Autophagy resulted from the formation of protein aggregates caused by AZA-dependent inhibition of DNMT2. Autophagy inhibition had an additive effect with AZA on AML cell proliferation and survival, increased ERE levels, increased pro-inflammatory responses, and generated immunogenic tumor-specific ERE-derived MAPs. Finally, autophagy was associated with a lower abundance of CD8+ T-cell markers in AML patients expressing high levels of EREs. This work demonstrates that AZA-induced EREs are degraded by autophagy and shows that inhibiting autophagy can improve the immune recognition of AML blasts in treated patients.

Breast cancer immunopeptidomes contain numerous shared tumor antigens.

Hormone receptor-positive breast cancer (HR+) is immunologically cold and has not benefited from advances in immunotherapy. In contrast, subsets of triple-negative breast cancer (TNBC) display high leukocytic infiltration and respond to checkpoint blockade. CD8+ T cells, the main effectors of anticancer responses, recognize MHC I-associated peptides (MAPs). Our work aimed to characterize the repertoire of MAPs presented by HR+ and TNBC tumors. Using mass spectrometry, we identified 57,094 unique MAPs in 26 primary breast cancer samples. MAP source genes highly overlapped between both subtypes. We identified 25 tumor-specific antigens (TSAs) mainly deriving from aberrantly expressed regions. TSAs were most frequently identified in TNBC samples and were more shared among The Cancer Genome Atlas (TCGA) database TNBC than HR+ samples. In the TNBC cohort, the predicted number of TSAs positively correlated with leukocytic infiltration and overall survival, supporting their immunogenicity in vivo. We detected 49 tumor-associated antigens (TAAs), some of which derived from cancer-associated fibroblasts. Functional expansion of specific T cell assays confirmed the in vitro immunogenicity of several TSAs and TAAs. Our study identified attractive targets for cancer immunotherapy in both breast cancer subtypes. The higher prevalence of TSAs in TNBC tumors provides a rationale for their responsiveness to checkpoint blockade.

BamQuery: a proteogenomic tool to explore the immunopeptidome and prioritize actionable tumor antigens.

MHC-I-associated peptides deriving from non-coding genomic regions and mutations can generate tumor-specific antigens, including neoantigens. Quantifying tumor-specific antigens' RNA expression in malignant and benign tissues is critical for discriminating actionable targets. We present BamQuery, a tool attributing an exhaustive RNA expression to MHC-I-associated peptides of any origin from bulk and single-cell RNA-sequencing data. We show that many cryptic and mutated tumor-specific antigens can derive from multiple discrete genomic regions, abundantly expressed in normal tissues. BamQuery can also be used to predict MHC-I-associated peptides immunogenicity and identify actionable tumor-specific antigens de novo.

Klf4 protects thymus integrity during late pregnancy.

Pregnancy causes abrupt thymic atrophy. This atrophy is characterized by a severe decrease in the number of all thymocyte subsets and qualitative (but not quantitative) changes in thymic epithelial cells (TECs). Pregnancy-related thymic involution is triggered by progesterone-induced functional changes affecting mainly cortical TECs (cTECs). Remarkably, this severe involution is rapidly corrected following parturition. We postulated that understanding the mechanisms of pregnancy-related thymic changes could provide novel insights into signaling pathways regulating TEC function. When we analyzed genes whose expression in TECs was modified during late pregnancy, we found a strong enrichment in genes bearing KLF4 transcription factor binding motifs. We, therefore, engineered a Psmb11-iCre : Klf4lox/lox mouse model to study the impact of TEC-specific Klf4 deletion in steady-state conditions and during late pregnancy. Under steady-state conditions, Klf4 deletion had a minimal effect on TEC subsets and did not affect thymic architecture. However, pregnancy-induced thymic involution was much more pronounced in pregnant females lacking Klf4 expression in TECs. These mice displayed a substantial ablation of TECs with a more pronounced loss of thymocytes. Transcriptomic and phenotypic analyses of Klf4 -/- TECs revealed that Klf4 maintains cTEC numbers by supporting cell survival and preventing epithelial-to-mesenchymal plasticity during late pregnancy. We conclude that Klf4 is essential for preserving TEC's integrity and mitigating thymic involution during late pregnancy.

Identification of antigenic epitopes recognized by tumor infiltrating lymphocytes in high grade serous ovarian cancer by multi-omics profiling of the auto-antigen repertoire.

Immunotherapeutic strategies aimed at enhancing tumor cell killing by tumor-specific T cells hold great potential for reducing tumor burden and prolonging survival of cancer patients. Although many potential tumor antigens have been described, identifying relevant targets when designing anti-cancer vaccines or targeted cell therapies remains a challenge. To identify novel, potentially immunogenic candidate tumor antigens, we performed integrated tumor transcriptomic, seromic, and proteomic analyses of high grade serous ovarian cancer (HGSC) patient tumor samples. We identified tumor neo-antigens and over-expressed antigens using whole exome and RNA sequencing and examined these in relation to patient-matched auto-antibody repertoires. Focusing on MHC class I epitopes recognized by CD8+ T cells, HLA-binding epitopes were identified or predicted from the highly expressed, mutated, or auto-antibody target antigen, or MHC-associated peptides (MAPs). Recognition of candidate antigenic peptides was assessed within the tumor-infiltrating T lymphocyte (TIL) population expanded from each patient. Known tumor-associated antigens (TAA) and cancer/testis antigens (CTA) were commonly found in the auto-antibody and MAP repertoires and CD8+ TILs recognizing epitopes from these antigens were detected, although neither expression level nor the presence of auto-antibodies correlated with TIL recognition. Auto-antibodies against tumor-mutated antigens were found in most patients, however, no TIL recognition of the highest predicted affinity neo-epitopes was detected. Using high expression level, auto-antibody recognition, and epitope prediction algorithms, we identified epitopes in 5 novel antigens (MOB1A, SOCS3, TUBB, PRKAR1A, CCDC6) recognized by HGSC patient TILs. Furthermore, selection of epitopes from the MAP repertoire identified 5 additional targets commonly recognized by multiple patient TILs. We find that the repertoire of TIL specificities includes recognition of highly expressed and immunogenic self-antigens that are processed and presented by tumors. These results indicate an ongoing autoimmune response against a range of self-antigens targeted by HGSC TILs.

OpenCustomDB: Integration of Unannotated Open Reading Frames and Genetic Variants to Generate More Comprehensive Customized Protein Databases.

Proteomic diversity in biological samples can be characterized by mass spectrometry (MS)-based proteomics using customized protein databases generated from sets of transcripts previously detected by RNA-seq. This diversity has only been increased by the recent discovery that many translated alternative open reading frames rest unannotated at unsuspected locations of mRNAs and ncRNAs. These novel protein products, termed alternative proteins, have been left out of all previous custom database generation tools. Consequently, genetic variations that impact alternative open reading frames and variant peptides from their translated proteins are not detectable with current computational workflows. To fill this gap, we present OpenCustomDB, a bioinformatics tool that uses sample-specific RNaseq data to identify genomic variants in canonical and alternative open reading frames, allowing for more than one coding region per transcript. In a test reanalysis of a cohort of 16 patients with acute myeloid leukemia, 5666 peptides from alternative proteins were detected, including 201 variant peptides. We also observed that a significant fraction of peptide-spectrum matches previously assigned to peptides from canonical proteins got better scores when reassigned to peptides from alternative proteins. Custom protein libraries that include sample-specific sequence variations of all possible open reading frames are promising contributions to the development of proteomics and precision medicine. The raw and processed proteomics data presented in this study can be found in PRIDE repository with accession number PXD029240.

Two types of human TCR differentially regulate reactivity to self and non-self antigens.

Based on analyses of TCR sequences from over 1,000 individuals, we report that the TCR repertoire is composed of two ontogenically and functionally distinct types of TCRs. Their production is regulated by variations in thymic output and terminal deoxynucleotidyl transferase (TDT) activity. Neonatal TCRs derived from TDT-negative progenitors persist throughout life, are highly shared among subjects, and are reported as disease-associated. Thus, 10%-30% of most frequent cord blood TCRs are associated with common pathogens and autoantigens. TDT-dependent TCRs present distinct structural features and are less shared among subjects. TDT-dependent TCRs are produced in maximal numbers during infancy when thymic output and TDT activity reach a summit, are more abundant in subjects with AIRE mutations, and seem to play a dominant role in graft-versus-host disease. Factors decreasing thymic output (age, male sex) negatively impact TCR diversity. Males compensate for their lower repertoire diversity via hyperexpansion of selected TCR clonotypes.

Induced pluripotent stem cells display a distinct set of MHC I-associated peptides shared by human cancers.

Previous reports showed that mouse vaccination with pluripotent stem cells (PSCs) induces durable anti-tumor immune responses via T cell recognition of some elusive oncofetal epitopes. We characterize the MHC I-associated peptide (MAP) repertoire of human induced PSCs (iPSCs) using proteogenomics. Our analyses reveal a set of 46 pluripotency-associated MAPs (paMAPs) absent from the transcriptome of normal tissues and adult stem cells but expressed in PSCs and multiple adult cancers. These paMAPs derive from coding and allegedly non-coding (48%) transcripts involved in pluripotency maintenance, and their expression in The Cancer Genome Atlas samples correlates with source gene hypomethylation and genomic aberrations common across cancer types. We find that several of these paMAPs were immunogenic. However, paMAP expression in tumors coincides with activation of pathways instrumental in immune evasion (WNT, TGF-ß, and CDK4/6). We propose that currently available inhibitors of these pathways could synergize with immune targeting of paMAPs for the treatment of poorly differentiated cancers.

Immunopeptidomic analyses of colorectal cancers with and without microsatellite instability.

Colorectal cancer is the second leading cause of cancer death worldwide, and the incidence of this disease is expected to increase as global socioeconomic changes occur. Immune checkpoint inhibition therapy is effective in treating a minority of colorectal cancer tumors; however, microsatellite stable tumors do not respond well to this treatment. Emerging cancer immunotherapeutic strategies aim to activate a cytotoxic T cell response against tumor-specific antigens, presented exclusively at the cell surface of cancer cells. These antigens are rare and are most effectively identified with a mass spectrometry-based approach, which allows the direct sampling and sequencing of these peptides. Although the few tumor-specific antigens identified to date are derived from coding regions of the genome, recent findings indicate that a large proportion of tumor-specific antigens originate from allegedly noncoding regions. Here, we employed a novel proteogenomic approach to identify tumor antigens in a collection of colorectal cancer-derived cell lines and biopsy samples consisting of matched tumor and normal adjacent tissue. The generation of personalized cancer databases paired with mass spectrometry analyses permitted the identification of more than 30,000 unique MHC I-associated peptides. We identified 19 tumor-specific antigens in both microsatellite stable and unstable tumors, over two-thirds of which were derived from noncoding regions. Many of these peptides were derived from source genes known to be involved in colorectal cancer progression, suggesting that antigens from these genes could have therapeutic potential in a wide range of tumors. These findings could benefit the development of T cell-based vaccines, in which T cells are primed against these antigens to target and eradicate tumors. Such a vaccine could be used in tandem with existing immune checkpoint inhibition therapies, to bridge the gap in treatment efficacy across subtypes of colorectal cancer with varying prognoses. Data are available via ProteomeXchange with identifier PXD028309.

Regulation of MHC Class I Expression in Lung Epithelial Cells during Inflammation.

Lung infections are a perennial leading cause of death worldwide. The lung epithelium comprises three main cell types: alveolar type I (AT1), alveolar type II (AT2), and bronchiolar cells. Constitutively, these three cell types express extremely low amounts of surface MHC class I (MHC I) molecules, that is, <1% of levels found on medullary thymic epithelial cells (ECs). We report that inhalation of the TLR4 ligand LPS upregulates cell surface MHC I by ∼25-fold on the three subtypes of mouse lung ECs. This upregulation is dependent on Nlrc5, Stat1, and Stat2 and caused by a concerted production of the three IFN families. It is nevertheless hampered, particularly in AT1 cells, by the limited expression of genes instrumental in the peptide loading of MHC I molecules. Genes involved in production and response to cytokines and chemokines were selectively induced in AT1 cells. However, discrete gene subsets were selectively downregulated in AT2 or bronchiolar cells following LPS inhalation. Genes downregulated in AT2 cells were linked to cell differentiation and cell proliferation, and those repressed in bronchiolar cells were primarily involved in cilium function. Our study shows a delicate balance between the expression of transcripts maintaining lung epithelium integrity and transcripts involved in Ag presentation in primary lung ECs.

CAMAP: Artificial neural networks unveil the role of codon arrangement in modulating MHC-I peptides presentation.

MHC-I associated peptides (MAPs) play a central role in the elimination of virus-infected and neoplastic cells by CD8 T cells. However, accurately predicting the MAP repertoire remains difficult, because only a fraction of the transcriptome generates MAPs. In this study, we investigated whether codon arrangement (usage and placement) regulates MAP biogenesis. We developed an artificial neural network called Codon Arrangement MAP Predictor (CAMAP), predicting MAP presentation solely from mRNA sequences flanking the MAP-coding codons (MCCs), while excluding the MCC per se. CAMAP predictions were significantly more accurate when using original codon sequences than shuffled codon sequences which reflect amino acid usage. Furthermore, predictions were independent of mRNA expression and MAP binding affinity to MHC-I molecules and applied to several cell types and species. Combining MAP ligand scores, transcript expression level and CAMAP scores was particularly useful to increase MAP prediction accuracy. Using an in vitro assay, we showed that varying the synonymous codons in the regions flanking the MCCs (without changing the amino acid sequence) resulted in significant modulation of MAP presentation at the cell surface. Taken together, our results demonstrate the role of codon arrangement in the regulation of MAP presentation and support integration of both translational and post-translational events in predictive algorithms to ameliorate modeling of the immunopeptidome.

PSMB11 regulates gene expression in cortical thymic epithelial cells.

The PSMB11 proteasomal subunit, expressed only in cortical thymic epithelial cells (cTECs), is essential for the development of functional CD8+ T cells. An attractive yet unproven theory holds that PSMB11 generates unique major histocompatibility complex class I (MHC I)-associated peptides required for positive selection. We recently reported that PSMB11 regulates the expression of hundreds of genes in cTECs, mainly by differential proteolysis of transcription factors. Thereby, PSMB11 maintains the distinctness of cTECs relative to medullary TECs (mTECs) and promotes cortex-to-medulla migration of developing thymocytes. These conclusions have been challenged by Ohigashi and colleagues, who suggest that their data show that PSMB11 uniquely controls antigen presentation without affecting cTEC biology. Here, we perform a comprehensive reanalysis of transcriptomic and proteomic data from the Ohigashi lab and confirm our original conclusions. This Matters Arising paper is in response to Ohigashi et al. (2019), published in Cell Reports. See also the response by Ohigashi and Takahama (2021), published in this issue of Cell Reports.

Beneficial autoimmunity improves cancer prognosis.

Many tumour antigens that do not arise from cancer cell-specific mutations are targets of humoral and cellular immunity despite their expression on non-malignant cells. Thus, in addition to the expected ability to detect mutations and stress-associated shifts in the immunoproteome and immunopeptidome (the sum of MHC class I-bound peptides) unique to malignant cells, the immune system also recognizes antigens expressed in non-malignant cells, which can result in autoimmune reactions against non-malignant cells from the tissue of origin. These autoimmune manifestations include, among others, vitiligo, thyroiditis and paraneoplastic syndromes, concurrent with melanoma, thyroid cancer and non-small-cell lung cancer, respectively. Importantly, despite the undesirable effects of these symptoms, such events can have prognostic value and correlate with favourable disease outcomes, suggesting 'beneficial autoimmunity'. Similarly, the occurrence of dermal and endocrine autoimmune adverse events in patients receiving immune-checkpoint inhibitors can have a positive predictive value for therapeutic outcomes. Neoplasias derived from stem cells deemed 'not essential' for survival (such as melanocytes, thyroid cells and most cells in sex-specific organs) have a particularly good prognosis, perhaps because the host can tolerate autoimmune reactions that destroy tumour cells at some cost to non-malignant tissues. In this Perspective, we discuss examples of spontaneous as well as therapy-induced autoimmunity that correlate with favourable disease outcomes and make a strong case in favour of this 'beneficial autoimmunity' being important not only in patients with advanced-stage disease but also in cancer immunosurveillance.

Atypical acute myeloid leukemia-specific transcripts generate shared and immunogenic MHC class-I-associated epitopes.

Acute myeloid leukemia (AML) has not benefited from innovative immunotherapies, mainly because of the lack of actionable immune targets. Using an original proteogenomic approach, we analyzed the major histocompatibility complex class I (MHC class I)-associated immunopeptidome of 19 primary AML samples and identified 58 tumor-specific antigens (TSAs). These TSAs bore no mutations and derived mainly (86%) from supposedly non-coding genomic regions. Two AML-specific aberrations were instrumental in the biogenesis of TSAs, intron retention, and epigenetic changes. Indeed, 48% of TSAs resulted from intron retention and translation, and their RNA expression correlated with mutations of epigenetic modifiers (e.g., DNMT3A). AML TSA-coding transcripts were highly shared among patients and were expressed in both blasts and leukemic stem cells. In AML patients, the predicted number of TSAs correlated with spontaneous expansion of cognate T cell receptor clonotypes, accumulation of activated cytotoxic T cells, immunoediting, and improved survival. These TSAs represent attractive targets for AML immunotherapy.

Most non-canonical proteins uniquely populate the proteome or immunopeptidome.

Combining RNA sequencing, ribosome profiling, and mass spectrometry, we elucidate the contribution of non-canonical translation to the proteome and major histocompatibility complex (MHC) class I immunopeptidome. Remarkably, of 14,498 proteins identified in three human B cell lymphomas, 2,503 are non-canonical proteins. Of these, 28% are novel isoforms and 72% are cryptic proteins encoded by ostensibly non-coding regions (60%) or frameshifted canonical genes (12%). Cryptic proteins are translated as efficiently as canonical proteins, have more predicted disordered residues and lower stability, and critically generate MHC-I peptides 5-fold more efficiently per translation event. Translating 5' "untranslated" regions hinders downstream translation of genes involved in transcription, translation, and antiviral responses. Novel protein isoforms show strong enrichment for signaling pathways deregulated in cancer. Only a small fraction of cryptic proteins detected in the proteome contribute to the MHC-I immunopeptidome, demonstrating the high preferential access of cryptic defective ribosomal products to the class I pathway.